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j3 ■ Abstract 
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Electron-positron pair production in an arbitrary polarized ultrastrong laser field is investigated 
in the first order perturbation approximation in which the Volkov states are used for convenient 
calculation of scattering amplitude and cross section. It is found surprisingly that the optimal 
pair production depends strongly on the polarization. For some cases of field parameters, the 
optimal field is elliptically polarized or evenly circularly polarized one, rather than the usual linear 
polarization as indicated by previous works. Some insights into pair generation are given and some 
interesting unexpected features are also discussed briefly. 
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I. INTRODUCTION 



Strong field physics has been an active research field since the invention of the technique 
chirped pulse amplification (CPA) in 1985 which draws the intensity of tabletop laser 
from gigwatt to erawatt and makes it possible to test the quantum electrodynamics 
(QED) theory in laboratories J3]. The Stanford Linac Acceleration (SLAC) experiment in 
1997 |3|, |4J have observed electrons-positrons (e + e~) pair production and revealed nonlinear 
QED effects, which have resurrected research interests in pair creation in strong fields. The 
forthcoming extreme-light-infrastructure (ELI) [5j will surely boost breakthrough in strong 
field physics research. This will be the strongest field in the world, though its intensity is 
approximately about 10 25 W/cm 2 and still 4 orders less than the Schwinger's critical field 
intensity 10 29 W/cm 2 [6j. 

Several methods in both of theory and numeric ca 



alyze pair production in external fields since 1930s jj-llG 



culation have been proposed to an- 



18| . One of the most important 



approaches is the perturbation method in which the corresponding Volkov states are used 



as 



13 



aases. It was pioneered by Volkov [7], and developed by Reiss (l9j, Nikishov 



20|, Midler 



16], Keitel 15|, |l6) and so on. But all of the above studies are based on particular polar- 



ized laser fields with either circular or linear polarization 
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20j. Most studies have 



manifested that the pair production rates are larger in linearly polarized fields compared 
than that in circular ones [bfl . ItI 18 . 21 ] in some conditions. Albeit these progresses which 



have been achieved, no complete analytical or numerical studies have been performed to 
reveal the influence of laser fields with arbitrary polarization on pair creation, especially in 
the framework of perturbation method even though this method seems simple and valid for 
the problem. The reason may be that relatively tedious derivation of scattering amplitude 
and cross section is involved in that case of elliptic polarization where the scattering partial 
waves have many coupled terms. 

Therefore, it is worthwhile to fill this research gap by exploring the relation between 
the elliptically polarized laser fields and pair production, which will give more insight into 
the SLAC experiment. Our study will show that a better pair production by elliptically 
polarized laser field can occur possibly. If we have got more information and recognition, it 
will be helpful not only to uncover more phenomena but also to regulate involved parameters 
in future experiments. 
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In the present paper, we will deal with one of these problems, through difficult theoretical 
and computing calculations. The case of head-on collision, occurring between a high-energy 
gamma photon and multiple photons of the laser fields, is considered. This process is called 
multi-photon reaction. Fortunately, we get the strict result in the framework of first order 
perturbation theory, and give thoroughly numerical calculations of some fixed parameters. 
There is, in fact, a surprising result that the pair creation relies strongly on the parameters 
of the laser fields, especially T (describing the elliptical polarization) in the condition that 



rjo (the Lorentz- invariant dimensionless strength parameter |22j) could be compared with t? 
(the collision energy). The optimal filed is not the linearly polarized field, but in some cases 
it is elliptically polarized and sometimes even the circularly polarized. 

The paper is organized as follows. Firstly, in Sec. II, we give a concrete theoretical deriva- 
tion of the total cross section in the first order perturbation theory using the Volkov states 
wave functions. In Sec. Ill, we present some demonstrations of numerical results which 
show clearly and vividly how pair creation is influenced by a few important parameters such 
as the polarization T, the Lorentz-invariant normalized field strength rjo, and the collision 
energy The main results will be summarized and some insights and discussions are given 
briefly in the final section. 

II. THEORETICAL FORMALISM 

In this section, we give a detailed derivation of the total pair production in the interaction 
between strong laser field and high ener gy p hoton. Some techniques used here are employed 
and developed from Ref. 22[ and Ref. 23J. First of all, some notations must be employed 
for simplicity and convenience. Throughout the paper, an elliptically polarized laser field, 
v4 M = a[ei cos {k ■ x) + Ts% sin {k ■ x)], is adopted, where a is the amplitude of vector potential 
A of laser field, k 2 = 0, k» = 0, 0, 1), e% = (0, 1, 0, 0), e% = (0, 0, 1, 0) and -1 < T < 1. 
The natural units, h — c—1, are also used. 

We define A = |(1 + T 2 ), the average number density of photons p(co) = 
(^[1 + (T 2 - 1) cos 2 (k ■ x))) = ^ and r, = ^ \{A^)\ = ^VX, where e and m 
are respectively the charge and the rest mass of an electron. 770 is defined as The scalar 
four-product is defined as A ■ B = A^B^ = AqB° — AiB % and^t = A^ , in which 7 M is the 
standard Dirac matrices. We also sometimes use <fi in place of k ■ x for brevity. 
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The solution of electrons/positrons under this type field could be well 
solved in terms of the known Volkov states [7, [22] , i.e. for electron it is 



(«) 



p.s 



1 + 



2{p-k) 



u(jp,s)exp[iSjpJk ■ x)] and for positron it is ^p-^ 



,(e+) 



Nfj 



2(p-k) 



v(p,—s)exp[iSfo2g\(k ' x )]- Here u(jp,s) and v(p,—s) are respectively 



the unit spinor of a free electron and a free positron, p and p the momenta, s and s the 



spins, and ST^ and S% are 



(e+) 
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(1) 



We will denote S(t} , simply as S^ e+ ^ and % as in the following discussion. Then, the 



(P.-s) 



(P. s ) 



amplitude for pair production process due to the laser field collision with a high energetic 
gamma photon of A' = Nqe' exp (— ik ■ x) is the following in the framework of first order 
perturbation theory: 



•r(«0 



XI 
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% = -ze / d 4 x 

= —ieNp-N^No J d 4 xu(pf, Sf)Mv(pi, — s$) exp iS^ e \x) — iS^ e \x) — ik' ■ x 
Through a detailed calculation, we have 
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= + - ^ M + eaQ • £i sin - TeaQ ■ e 2 cos 
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Here, we have defined the following quantities: 
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If 
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_2i_ 
qi-k" 1 



-z sin 
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Thus, we find the amplitude 

S fi = -ieN Pf N Pi N f d A x u(pf,Sf)Mv(pi,-Si)x 

exp {i [q f;ii + g i;jU - k^] - iz sin (0 - <f> ) - iy sin (20) } 

oo oo 

= -ieN Pf N Pi N J m{y) E (2vr) 4 5 (g /ijU + - ^ - 2mk fl - nk^) 

m=—oo n=—oo 

u(p f ,8 f ) [B n (z)M + B ln (z)M l + B 2n (z)M 2 + B 3n (z)M 3 ] v(p h 
where some notations have been adopted as 



M = M + Mi cos + M 2 sin + M 3 cos (20), 
M = £ - Ae 2 a 2 (fc ' £ ^ 

< Mi = ea 



M, = Tea 



2(p r k)(p i: k)- 
2(p/-fe) 2(pi-fc) 



2(p / -fe) 2(pi-fc) 

M 3 = (A - l)e 2 a 2 (fcV) * 



and some special functions like J n (y) and B n (z) etc are given by 
exp [-iy sin (20)] = X Jn(v) ex P H2n0], 

?1= — oo 

oo 

exp [-iz sin (0 - O ))] = XI ^n(^) exp [-incp], 

n=— oo 

oo 

cos 0exp [— sin (0 — O )] = ^ B ln {z) exp [— m0], 

n=— oo 

oo 

sin exp [— sin (0 — O )] = XI B 2n (z) exp [— m0], 

n=— oo 

oo 

cos (20) exp [— zzsin (0 — O )] = X B 3n {z) exp [— in0]. 

n=— oo 

The amplitude can be also rewritten in a more compact form, 
S fi = -ieN Pf N p .N x 

oo oo 

X Jm(y) X (2vr) 4 5 (g /;M + - ^ - 2mk tl - n/c^) Sf)M n v(p u -Si 



(3) 



(4) 



(5) 



where 



M n = B n (z)M + B ln (z)M 1 + B 2n {z)M 2 + B 3n {z)M 3 . 



What we are interested in is |S/j| 2 . With a deliberate consideration on Eq. (j5]), the main 
task has been converted to how to tackle Mfi >n = u(pf, s f)M n v(pi, — s«). Then we will have 
to deal with Mf^ n M% m in the general elliptic polarized field rather than Mfi >n M% n in the 



circular polarized field 22] because there are two series summations existed in Eq. fl5]) in 
the case of general elliptic polarization. 

The final total cross section a is obtained by averaging over the photon polarization A 



and summing over the e+e spins, M fijn M* i m = \ J2 M fi,n M }i 

A Si,Sf 



M fi7n M 
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\ E Tr 
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hT.Tr [(fa + m )M n (p i - m^M^ ] . 
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Before proceeding, we define some notations as below 



e n = B ln {z)e 1 + TB 2n (z)e 2 , 
e n = B* ln (z)e 1 + TB* 2n (z)e 2 , 



(7) 



and 



" A{q r k){qi-k) ' 

s = 2(k' ■ k), 
m l — m o + e 2 a 2 A. 



(8) 



Then we can easily get k ■ e n — 0, since k ■ E\ — and • e 2 = 0. 



We should deal with Mf^M^ and Mf itm Mj i<n together, based on that the final 
result of \Sfi\ 2 should be real, thus we can define a convenient quantity M| i;mn = 
\ (Mfi tn M% m + Mfi iTn M% n ) . After a tedious and careful calculation, we finally obtain 

jh n \j2Tr [{P f + m )M n (P l - m )7° 'Mfrf] 

(9) 



Ml 

ji;mn 

" K. A 

+ £ Tr [(^ + m )M m (^ - m ) 7 °Mt 7 ] 

A 

E m n,0 E mn ^\ -\- E mn2 E mn ^ ) 

where E mnfi , E mn>l , E mn>2 and E mnt3 are defined as 

^mn,o = coa[(n-m)<j> ]J n (z)J m (z) 

+rf cos [(n - m)0 o ] U(^) J m (^) - Jm+l(2) /" +l(2) - (1 - 2m), 



E m n,l 



4m 2 , 



(2N — m — n) cos [(n — m)0 o ] J n (z)J m (z), 



E„ 



mn,2 



[1 - \)rj 2 cos[(n -m)0 o ] cos (2(f) )J n (z)J m (z) x 



[l-2u] [^-f m|J m (*)|f ln|J n (z 



n 2 +m 2 



l_d_ 

z dz 



In \JJz) J m (z 



and 

E mn ,3 = -(1 - \)r] 2 sm[(n - m)(f)o]sm(2(f) )J n (z)J m (z)x 

{(n-m) [±MJ n {z)J m {z)\-^]-^^}. 

It is necessary to use the basic Bessel relations and the equation, pfpi = 7;Ns — 2i] 2 m 2 l u—m 2 ) , 
to obtain Eq. 

Finally, we get the differential cross section, 
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mi,ni=-oo m2,ri2=— oo 

(10) 

where p w is defined as ^f^, = 1 when a = b and when a ^ b. 

To derive the total cross section a of pair production, we will first only consider the 
interaction between the gamma photon and photons taken from the laser field, in which 
If + Qi = k' + Nk is satisfied. To deal with Eq. ( ITOj) . we work in the center coordinate of 
the pair created. Then the following quantities could be represented as 
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q 

CO 



(Qf,Qf) = (q°,q), 
(q°,Qi) = (q°,-q) 

ml + A 2 e 2 a 2 
Nco, q = co . 



(11) 



and 



s 

q r k 
u 

cos 9 



Auo co = 4Nu 2 



to oo 



\q\co cos 9, 



Ns 



(12) 



4 <jg — | q\ 2 cos 2 



Here, we first write a useful integral, J J d 3 qfd 3 qi-^)5 + — k„ — Nk^j, which in 
fact could be treated as an integral operator, into a more compact form. With the help of 
Eq. (j2J), Eq. ([Til . Eq. ( fl"2l and the property of Dirac function 5, we have 

/ / d 3 q f d 3 qi^5 (q f . !fl + q i]fl -k'^- Nk^) = f d 3 q^5 (2q - to' - Nio) 



J dip J dcos9 J d \q\ ^-5 (2q - oo' - Nco) = j dip j dcos9 

h [ dip [ du ; 1 

2 J ^ J uy/u(u-l) 



2?o 



Then the total cross section <jn of the pair production, in the interaction between one high 
energy gamma and iV photons of the laser field, can be presented in a compact form, 

2a \ C dip C N du / x 



a N 



E 

mi,m=— oo m2,ri2=— oo 



Jm 1 (y) J m 2 (y)] , 



where a is the fine-structure constant, $ 



and 



4(1+J? 2 ) 



with N > n because the 



minimum photon number no is need to overcome the energy gap to create the e + e pair. 



+00 



The total cross section a is Yl a N- It is easy to get from Eq. fTTTl) and Eq. (Tl2]) the 



N=n 



relation 



iV.s 



< U < 



-, i.e. u G 1, u^] 



It is noted that Eq. (fl3|) is determined by only three parameters, T, rj and 1?. The 
remaining problem is how to handle </>o, ?/ and 2. It is possible to get their final expressions 
from the relations in Eq. (J2J), Eq. ([8]), Eq. (ITT]) and Eq. (fl2|) . After a careful calculation 
with these relations, we finally obtain 
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As a check we get easily the result for the circular polarization, 
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where no is the minimum integer among possible values for N that satisfy the condition 
Un > 1.0. It is also the minimum number to overcome the energy gap for pair production 
mentioned above. Obviously our results for this special case is the same as that in Ref. 22 ]. 



III. NUMERICAL RESULT 



Based on the theoretical formula obtained in last section we can calculate numerically the 
cross section for the pair production for different parameters. The C code of Bessel function 



in Ref. 



24j is used in our program. 



With Eq. (|T3|) . the corresponding cross section in MKSA units should be 

/ \ 2 „ 2tt u n 

a N = ( JL) f dtp du / x 

oo oo _ \ / 

S S ^n 1 +2mi,Af^n2+2m 2 ,iV [^/j ;ni „ 2 ^mj (?/) ^m 2 (Z/)] , 

mi,ni=-oo m2,fi2=— oo 

where N >tiq. 

Our numerical simulations are based on Eq. f|T5|) . and the numerical results have been 

2 fi 2 + °° 

scaled by ° 2 2 , then the total cross section a is ^2 &n- The main numerical results will 

m ° C 7V=1 

be demonstrated graphically in the following subsections. It is noted that we have also 
reproduced the numerical result of Ref. {22] for the special case, i.e. the circular polarization. 



A. Pair production cross section dependence on laser field polarization and inten- 
sity 

We consider first the laser field applied in the SLAC experiment and give an insight 
about the experiments with d = 1.0511647 and r] = ea/m = 1.0. Note that it is different 
from r] = r/ov^A in our general case since the laser field polarization affects A. The relation 
between the cross sections of pair production and the parameter T is displayed in Fig. [TJ 

Surprisingly we can clearly see that the optimal polarization parameter T is 1.0 or —1.0 
rather than as thought by some people. On the other hand, the total production cross 
sections a at T = 0.3, 0.4 and 0.5 are also very close to the optimal one. Furthermore, in the 
case of linear polarization the total pair production cross section is the smallest. It can be 
easily obtained that a(T = l)/o"(T = 0) rj 1.34, and a(T = 0.3)/a(T = 0) » 1.31. These 
ratios indicate that the parameter T plays a very important role in directing the experiments 
to get more obvious and easy-detected results. We can also see that the numerical result is 
symmetric with respect to T = as is expected by the theoretical requirement in Eq. ffl~5l) . 

Then we will unearth how the cross sections, and a, are influenced by 770, the intensity 
of the polarized laser field. Usually the parameter $ is kept as a constant in experiments 
when rjo is changed. Now we will calculate the pair production cross sections for different 
770. Here the strong constraint, un = ^1^2^ > 1, should be taken into account seriously, 
otherwise the integral of Eq. (TT5T) would be zero. Given •& = 1.0511647 some numerical 
results with r] = 0.1, 0.5 and 1.5 are plotted in Fig. El Fig. [3] and Fig. HJ respectively. We 
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-1.0 -0.5 0.0 0.5 1.0 

Y 

FIG. 1: (Color online) The pair cross section <jn for different number of photons dressed, i.e, from 6 
to 30 and a, the summation of them, are plotted. Here, for example, N = 6, represents the number 
of photons dressed from the laser field in the multi-photons process, rjo and i? are the parameters 
fixed in the numerical simulation, while T is controlled. We only plot for N € {6, 7, . . . , 15}, and 
a « ^jv=6 because an = when iV = {1, . . . , 5}. All other figures will have these adoptions 
unless otherwise specified. 

also give the total cross sections a for a series of different field intensities in Fig. [5] in order 
to see the effect of r) . 



TABLE I: The relation between the optimal T which makes maximum and N when 770 = 1.0 
and 1? = 1.0511647. "0.9/1.0" means that the cross section ajsr when T = 0.9 and T = 1.0 cannot 
be distinguished in our numerical simulation. Similar in Tables IITI and iDIl 



N 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


|T| 


0.1 


0.3 


0.4 


0.5 


0.6 


0.7 


0.9/1.0 


1.0 


1.0 


1.0 



As an example, in Table U we list the values of |T| corresponding to the optimal pair 
production cross section ctn for different N when d = 1.0511647 and 770 = 1.0. Actually, we 
have studied all the data we have got from numerical calculations, and find a fact that there 
exists a tendency: the larger N, the greater |T|. This means that the polarized field becomes 
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FIG. 2: (Color online) Pair production ajy and a when r]Q = 0.1 and i9 = 1.0511647. Here 
a « Sjv=4 and for N £ {7,... ,10} have been omitted for their insignificance. 

closer to the circular one for the optimal pair production as N increases. This behavior 
is more obvious for smaller r] where N is larger. Therefore, a dominant contribution to 
arises from |T| = 1 when iV is large. However, this conclusion for total pair production is not 
valid in general because the iV-order cross section attributed to the total cross section 
decreases as N increases. This point is also seen in Fig. [5j for example, where the optimal 
total pair production for t]q = 1.5 is not at |T| = 1 but at 0.4. It should be emphasized that 
for the total cross section there is a trade-off between i] and T which makes a optimal. 

TABLE II: The relation between the optimal T which makes a maximum and 770 when 1? = 



1.0511647. 



V0 


0.1 


0.3 


0.5 


0.7 


1.0 


1.2 


1.5 


|T| 


0.1 


0.1 


1.0 


1.0 


1.0 


0.3/0.4/0.5 


0.4 



On one hand, for smaller 770 the main contribution to total a comes from a no since 
other terms (o"Ar >no ) decrease rapidly, see Fig. [2j In this situation the optimal polarization 
parameter approaches the linear one, see also Table IITI in the cases of 770 = 0.1 and 0.3. On 
the other hand, as ^0 increases to ~ 1.0, the field polarization parameter of optimal total 



1 ' 1 1 1 1 1 1 r 

r, =0.1,9=1.0511647 
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TABLE III: The relation between the optimal T which makes a maximum and $ when t]q = 1.0. 



r) 


1.0511647 


3.0 6.0 10.0 


15.0 21.0 28.0 36.0 


45.0 


55.0 


|T| 


1.0 


0.5/0.6 0.6 0.3 


0.3 0.3 0.3 0.2/0.3 


0.2/0.3 


0.2/0.3 



pair production cross section becomes the circular one, |T| = 1. The same results are shown 
in Tabled for very large N, because in these cases the terms of &N >no decrease more slowly 
so that the whole contribution from these large N terms exceeds that from those small N 
terms. This point can also be seen clearly in Fig. [1] and Fig. |3j However, as rjo increases 
further a typical nonlinear feature arises between |T| and t]q for the optimal pair creation. 
From Fig. [4] and Table ILTI the trade-off between |T| and t]q contribution to total cross section 
compels the optimal polarization to locate at about |T| = 0.4, which is between linear and 
circular polarizations for our studied field parameters. 

In a word, although the pair production cross section increases monotonically with laser 
field intensity, its optimal value depends strongly on laser field polarization, see Fig. [5j The 
former can be understood by the perturbation concept but the later can be recognized as 
the typical nonlinear interaction between fields and vacuum in QED, i.e. the laser dressed 
electrons/positrons colliding with a high energy photon. 

B. Pair production cross section dependence on i? 

Besides two important parameters, T and 770, the pair creation is also closely related to 
the parameter The problem, in fact, is how to tune the parameter to improve the pair 
production, given that un = 47^2) > 1- The idea is conformed with the physical picture as 
V$ = a/2(&' • fe)/mo, which means the collision energy of two photons. Usually one believes 
that the larger $ is, the more likely the process of pair generation occurs. But things may 
not be so simple. We will find that the pair production will surely be greatly enhanced, in 
the first order, then begin to decrease in quantity when d is bigger than a critical i? c , which 
again exhibits the nonlinear characteristic of QED process. Some values of •& are selected 
for fixed r] , which are 3.0, 6.0, 10.0, 15.0, 21.0, 28.0, 36.0, 45.0 and 55.0 and the numerical 
results with d = 10.0, 15.0, 55.0 are demonstrated graphically from Fig. [6] to Fig. [HJ We 
also give all the total cross sections a for different $ in Fig. [9] and Fig. [10] in order to make 
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FIG. 3: (Color online) Pair production cr^ and cr when t]q = 0.5 and i9 = 1.0511647. Here 
a X^]v=5 and ajq for N € {5, 11, ... , 14} have been omitted for their insignificance. 
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FIG. 4: (Color online) Pair production and a when tjq = 1.5 and i? = 1.0511647. Here 
a « 5Zat=9 ^JV and cj at for N G {19, . . . , 40} have been omitted, but the contribution from them is 
not insignificant. 
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FIG. 5: (Color online) All the total cross sections a for different ryo when •& = 1.0511647 are plotted 
in one figure to unearth more information. 

convenient comparisons with each other. 

Now we will see how d influences the pair creation. Firstly, from Fig. |9j we can see that 
the pair production is greatly enhanced with the increase of $ for fixed ?7o, same behavior as 



shown in Ref. [22j for the circularly polarized field. But it is not monotonically increasing 
with d since the total cross section o starts to decrease when $ > 15.0 in our numerical 
results. The physical mechanism of this phenomenon seems difficult to grasp, however, a 
closer inspection of the integral of Eq. ( fTBI) shows that there exist two mathematic reasons 
for this counterintuitive phenomenon, i.e. the influence of Bessel functions and the inverse 
relation between the total cross section a and $ in the coefficient ^ . Certainly whether this 
phenomenon is kept in multiphoton perturbation regime, in which the higher order terms are 
included, or even in nonperturbation regime, is still an open problem. This difficult problem 
may be overcome in future possible theoretical calculations as well as more experiments. 

a is a summation of different when N is subject to > 1.0. Let us remember a 
fact that if the main contribution comes from the first term with least number N = n then 
the optimal pair production is prone to the linearly polarized field. This is more obvious for 
larger as shown in Table IIHI But there are some exceptions for smaller for example, 
|T| = 1 corresponds to parameter $ w 1.05 in SLAC experiment case. More interesting, in 
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FIG. 6: (Color online) Pair production ct/v and a when tjq = 1.0 and $ = 10.0. Here a sa J2n=i a N 
and cT/v for N 6 {6, . . . , 10} have been omitted for their insignificance. 
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FIG. 7: (Color online) Pair production ct/v and a when r/o = 1.0 and i? = 15.0. Here a ~ X^iv=i a N 
and CTjv for G {6, . . . , 10} have been omitted for their insignificance. 
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FIG. 8: (Color online) Pair production ct/v and a when tjq = 1.0 and $ = 55.0. Here a sa J2n=i a N 
and aj\r for N G {5, . . . , 10} have been omitted for their insignificance. 

our first order perturbation theory, when i? is large the pair production seems insensitive 
to the polarization parameter T. From our numerical results we can see clearly, in Fig. 
[U Fig. [9] and Fig. [TDJ that on one hand there is relatively large fluctuation of the cross 
section with respect to the polarization parameter for •§ £ [1.05, 21.0] and on the other hand 
there is indeed relatively small change of the cross section with respect to the polarization 
parameter T when d > 28.0. Although the main reason is unclear physically, the conclusion 
from numerical results is still very meaningful and useful. It provides an intuitive conclusion 
that when the collision energy of two photons is large enough usually one can directly choose 
the linearly or circularly polarized laser field, which is relatively easier to get in laboratories, 
for the pair production. With such special choices of polarization the cross section of the 
pair production does not differ from the optimal one greatly (less than 15% when •& > 36.0, 
as shown in Fig. [TUl) . 

IV. DISCUSSIONS AND CONCLUSIONS 

We have obtained the rigorous formula in the first order perturbation theory and per- 
formed thorough numerical computations for the pair production. Our numerical results for 
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FIG. 9: (Color online) All the total cross sections a in different i? when rjQ = 1.0 are plotted in one 
figure to unearth more information. Here i? = 1.0511647 is shortened as i? = 1.05 in the figure. 




FIG. 10: (Color online) This figure is corresponding to Fig. [9] and here "Ratio" is the ratio of ay 
and the corresponding a of the circular in the same parameters $ and tjq. We have omitted the 
cases when i? = 1.0511647, 3.0. 
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the parameters of the SLAC experiment have clearly manifested that the linear laser field 
is not the optimal one, but rather the circular one, although the partial cross section <r 6 in 
the case of linear polarization is larger than that in the case of circular one. Surprisingly, 
neither linear nor circular but the elliptically polarized field corresponds to the optimal pair 
creation cross section in some cases shown in Table [III and Table IIHI 

Our numerical results have also demonstrated that the leading contribution in the pair 
production comes from the first term subject to un = 4(i^ 2 ) > 1-0, when the polarization is 
nearly linear; moreover, things begin to change as larger N is considered, since ex at for larger 
N is prone to the circular one. It might be the main reason that causes the balance between 
the linearly and the circularly polarized and also the reason of many nonlinear dependence 
of pair production on other system parameters. From our numerical results, we also see 
that when 770 is large the contributions from <jn with large N are very important. While 1? 
is large the contribution from with large N are not so important. For example, given 
1] = 1.0, when we calculate the total cross section for $ = 1.0511647, at least Yl%=e^N is 
used; however, for $ > 6.0, summation by ^2n =1 <Jn as an approximation of a is already 
good enough. By the way for i] = 1.5 and $ = 1.0511647 a calculation by Y1n=9°'n is 
needed. 

We have exhibited the abundance phenomena in the process of photon-multiphoton re- 
action that the laser field polarization plays a key role in optimal pair production when the 
other parameters are fixed, which can be seen from Fig. [T] to Fig. [10] and three tables. 
Our work may provide some understanding about some previous experimental phenomena 
{3I, [2l| and also direction for future experiments. Although a great effort has been made 
to understand the polarization effect on the pair production, there still remain some open 
problems to be solved, especially the physical mechanism behind the complex calculational 
formulae. We conjecture that there may be nonlinear resonance or/and interference terms 
appeared in different modes of scattering matrix elements for the transition amplitude which 
would enhance or reduce the pair production cross section in different polarization cases. 
However, more theoretical and experimental research are still needed in the future to have a 
deeper understanding of vacuum decay in ultrastrong laser field to create electron-positron 
pairs, which serves as an important test for nonlinear QED. 
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